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Abstract 

The only source of ice formation in the Arctic during summer is a layer of ice be- 
tween an under-ice melt pond and the underlying ocean, called false-bottoms. The 
simultaneous growth and ablation of false-bottoms is governed by both of heat fluxes 
and salt fluxes. This is a two-phase Stefan problem with two free boundaries. We first 
use Green functions to reduce this problem to solving a system of nonlinear integral 
equations, and then apply the contraction principle to prove the existence and unique- 



(~| , ness of the solution under suitable data. 
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1 Introduction 

>, 

^^ ■ There exists fresh water in the Arctic during summer. First, melt water collects in 

^s^ , surface melt pond (melting under the sun) which is the most important reservoir. Second, 

CN ' melt water can percolate into the ice matrix to form an under-ice melt pond. At the interface 

t~^ ■ between this fresh water and the underlying salt water, double-diffusive convection of heat 

^^ , and salt occurs, leading to the formation of underwater ice called "false-bottoms" (see, e.g, 

Q I [Hlllll])- Note that salt water has the double properties: 

i) it does not freeze even for temperature < 0*^0, 

ii) it dissolves ice when it is in contact with ice. 



Surface melt pond 
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Sea water (Oc) Under ice melt pond (Wa) Ice forming (False -bottoms) 

Figure 1: Ice formation in the Arctic during summer. 

Thus false-bottoms is a layer of ice, which relatively protects the fresh water from being 
mixed into the underlying salt water. Such false-bottoms is the only significant source of 



ice formation in the Arctic during summer. The main difference between the false-bottoms 
case and the normal sea ice case is the presence of the growth at the upper interface of the 

first one. 

Now we consider a one-dimensional model describing the simultaneous growth and 

ablation of the ice of false-bottoms. Here we have three environments: the ocean (Oc), the 

ice of false-bottoms (Fb) and the fresh water (Wa). Denote by T{x,t), S{x,t) the 

temperature and the salinity, and denote by hQ{t), hu{t) the free boundaries at the 

interfaces ice-ocean (Fb-Oc) and ice- water (Fb-Wa), respectively. 




Figure 2: The one-dimensional model. 
At the interface Fb-Oc, we apply the first principle of thermodynamics, i.e., 

(variation of energy = variation of heat flux through the interface) . 
The net amount of heat transferred through the interface Fb-Oc in a section s is equal to 

At/ = —dhop^LfS, 

where pi is the density of the ice and Lj is the latent heat of fusion. On the other hand, 
the difference of the heat fluxes through a section s in the ice and the ocean during a time 

dt is 



A$ 




Fb,ho 



, dT 
ox 



sdt, 



Oc,ho. 



where A/, Xq are thermal conductivities of the ice and the ocean. Finally, the law of 
conservation of energy, i.e. AC/ = A$, leads to the Stefan condition for the heat balance 

at the interface 



h'oit) = \i 



dT 
dx 



Xo 



FbM) 

where 



dT 
dx 
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A/ = — -f—, Xo 
PiLf 
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For simplicity, we can neglect the salt of the ice of false-bottoms. The water near the 
interface Fb-Oc is a mixture of melt water, which melts from the ice of false-bottoms, and 

sea water. This water freshens at the rate 

So{t)h'o{t), 

while salt diffuses into this water at the rate 
dS' 



-D 



dx 



Oc,h(i 



where 5'o(t) = S{hQ{t),t) is the salinity of the ocean at the interface and D is the 
molecular diffusivity of salt in sea water. The balance of salt at this interface leads to the 

conservation condition 



Soit)K{t) = -D-^ 



(1.2) 



Oc,ho 



The establishment of these equations bases on Martin and Kauffman |4j and Notz et al. 

m. 

At the interface Fb-Wa, we use the simplified scheme in [3], in which the temperature of 
the water in the under-ice melt pond is kept at 0°C. In particular, it leads to the 

boundary condition at the interface 



T{hu{t),t)=0. 



:i.3) 



Furthermore, due to the neglect of heat fluxes from the fresh water above the 
false-bottoms, a thermodynamic condition similar to (jl.ip reduces to the Stefan condition 

at the upper surface 



Kit) = Xj 



dT 

dx 



:i.4) 



Fb,hn 



Equations (ll.ip . (ll.2p . (ll.4p form a closed system in connection with the diffusion 
equations for heat and salt in the ice and in the ocean 
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, hu{t) > x> ho{t) 

Fb 


d'^T 
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, ho{t) > X > — cx) 

- Oc 


and 


5251 

dx^\ 


, ho{t) > X > — CX), 

Oc 



(1.5) 
(1.6) 

(1.7) 



where Dj and Dq are the thermal diffusivity. Here all constants are positive. 

The problem is of finding {ho, hu,T, S) from the system (ll.ip - (ll.7|) . where T = T(x, t) 

with hu{t) > X > — CX) and S = S{x,t) with /io(i) > x > —ex. 



2 Existence and Uniqueness 

We want to give the existence and uniqueness of the solution {ho, hu,T, S) under 
suitable data. For this purpose, it is reasonable to require the initial conditions /io(0), 
hu{0), T{x,0) = T°{x) with /i„(0) >x> -oo, and S{x,0) = S^{x) with ho{0) > x> -oo. 
Moreover, assume that the freezing temperature T(/io(t),t) = 7o(^) with < t < o" has 

been known . 
For rigour, we need some assumption on the data as well as a definition of the solution. 

Hypothesis 2.1. Let a > 0. Assume that 

(HI) ho{0) < hu{0). 

(H2) T^{x) is continuous at x = hQ{0) , x = hu{0) — , and satisfies the compatible conditions 

T^{ho{0)) = To(0), T^{hu{0)) = 0; T^(x) is continuous and bounded in (— cxd, /io(0)) and 

{ho{o),hum. 

(H3) S^{x) is continuous and bounded in (— cxd, /io(0)). 
(H4-) TQ{t) is continuously differentiable in < t < a. 

Definition 2.1. We say (ho, /i„, T, S) is a solution in < t < a if 

(CI) ho^t), hu{t) is continuously differentiable and hu{t) > ho{t) in < t < a. 

(C2) Tt, Txx is continuous in < t < a, hu{t) > x > ho{t) and ho(t) > x > — oo; T(x,.) is 

continuous att = 0; T{.,t) is continuous at x = hQ{t) and x = hu{t) — ; Tx{-,t) is continuous 

at X = /io(t)— , X = /io(t)+ and x = hu{t) — . 

(C3) St, Sxx is continuous in < t < a, ho{t) > x > — oo; S(x,.) is continuous at t = 0; 

S{.,t) is continuous at x = ho^t). 

(C4) The system lll.l\) - [T7}\ ) is satisfied inO < t < a with the initial conditions (/io(0), hu{0),T^{x), S^{x)) 

and the freezing temperature T(/io(t),t) = To(t). 

Our main result is as follows. 

Theorem 2.1. Assume that (/io(0), /i„(0), r°(x), S°{x),To{t)) satisfy Hypothesis (H1)-(H4) 
for some cr > 0. Then the system /il.l]) - p7?\ ) has a unique solution {ho,hu,T, S) in < 
t < a* for some a* G (0,o"]. Moreover, this solution can be extended uniquely whenever the 
condition hu{a*) > ho{cr*) still holds. 

Our problem is a two-phase Stefan problem with two free boundaries. We shall follow the 

approach of A.Friedman ([2], Chap. 8), which deals with a classical one-phase Stefan 

problem with one free boundary. We first reduce the problem to solving a system of 

integral equations with respect to So{t) = S{ho{t) — ,t), vi{t) = Tx{ho{t) — ,t), 

V2{t) = Tx{ho{t)+,t) and V3{t) = Tx{hu{t) — ,t), and then solve this system by contraction 

principle. 

2.1 Preliminaries 

Let a{t),b{t) be continuously differentiable functions and b{t) > a{t), t>0. Let k > 
be a constant and let u{x, t) be a solution of the heat equation 

du d u 

with the initial condition m(x,0) bounded in (a(0), ^(0)). 
Introduce the Green's function of equation (j2.ip . 

G{x, t; C, r) = — ^ ' exp ' ^ ^^ 



2^7rK(t - r) V Mt - -r] 
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where H is Heaviside function, 



H{t) 



1, if t >0, 
0, if t < 0. 



We have the following Lemma to represent the solution of (j2.1|) from its initial condition 

and boundary conditions. 

Lemma 2.1. For t > and a{t) < x < b{t) we have 
t 

U{x,t) = G{x,t;b{T),T)[KU^{b{T)-,T)+u{b{T),T)b'{T)]dT 



t 
- / G{x,t;a{T),T) [Ku^{a{T)+,T) +u{a{T),T)a'{T)]dT 



K / G^{x,t;b{T),T)u{b{T),T)dT + K / G^{x,t;a{T),T)u{a{T),T)dT 





m 



+ I G{x,t-^,0)u{C,0)dC. 



a(0) 



Proof. Note that 



Gr + kG^^ = for all t < t, and G{x, t; S,, t—) = 5{x — ^), 

where 5 = H' is Dirac delta function. 
Integrating the Green's identity, here u = u{(^,t), 

over the domain a{T) < S, < 5(t),0 < r < t, we will obtain the desired result because 

t b{T) 



(2.2) 



d ( du dK\ 

a{T) 




9e " " 9e . 



dT 



C=a(r) 



t t 

G{x,t;b{T),T)u^{b{T)—,T)dT— / G{x,t;a{T),T)u^{a{T)+,T)dT 











i t 

G^{x,t;b{T),T)u{b{T),T)dT + / Gg(a;, t; a(r),r)u(a(r),r)dT, 





and 



t Kr) 



J J ^{uG)dCdr 



a(T) 



d_ 



b{r) 



b{r) 



U 
V(r) 

T=t- 



GdS, - [uG]^=b(^)6'(r) + [uG]^=a{T)a' (r) 



> dr 



uGdi 



a{r) 



T=0 



t t 

wG']^=b(^)6'(r)dr+ / [uG]^=„(^)a'(r)dr 



t 



m t 

u{x,t)- / G{x,t;^,0)u{^,0)d^- j G{x,t;b{T),T)u{b{T),T)b' {T)d7 



a(0) 

t 

+ / G{x,t;a{T),T)u{a{T),T)a'{T)dT. 


Here we have made use the following equality, with v = uG, 

b{r) I Hr) 



j §^i<^'^))dC = ^ J v{C,T)dn -v{h{T),T)h'{T)+v{a{T),T)a'{T). 



air) \.{t) 

In fact, in the case a{t) = we have 

6(r) 



J ^ {v{^, T))d^ = j H{b{T) -0^ {vi^, r)) d^ 



— [Hibir) - Ov{t r)] dC- g:^ mKr) - 0] <?, r)dC 



(oo \ oo 

I H{b{T) - i)v{i, T)di j - j 6{b{T) - Ob'{r)v{^, T)d^ 



d_ 
dr 



.0 



V 



vi^,T)dn -v{biT),T)b'iT). 



If a{t) is not constant, we can write 

b(r) b{r) 



6(r) 6(r) a{T) 

I |-(<e,r))d^= I ^ivi^,T))d^- J -^ivi^,T))d^, 



a{T) 

and the desired equality follows. This completes the proof. 



D 



Remark 2.1. The result in Lemma \2.1\ still holds for b{t) = +00 or a{t) = —oo. For 
example, if a{t) = —00 then the fomular in Lemma \2.1\ reduces to 

t 

U{x,t) = G{x,t;b{T),T)[KU^{b{T)-,T)+u{b{T),T)b'{T)]dT 



t 6(0) 

-K I G^{x, t; b{T),T)u{b{T),T)dT + I G{x, t; ^, 0)u(^, 0)d^. 

-00 

We shall need the following Lemma (see [2], page 217, Lemma 1). 

Lemma 2.2. Let p{t) be continuous and let s{t) > satisfy the Lipschitz condition, < 
t < a. Then, for < t < a, 

t t 

lim K / p{t)Gx{x, t; s{T),T)dT = -p{t) + k p{t)Gx{s{t),t; s{t), T)dT. 

x->s{t)- J I J 



Remark 2.2. This lemma gives the jump relation at the boundary s{t). In applications 
later, sometimes we need to note that Gg = —G^- Moreover, for the right limit we have 

t t 

lim K p{t)Gx{x,t;s{T),T)dT = --p{t) + K p{t)Gx{s{t),t;s{T),T)dT. 

x-*s(t)+ J 2 J 



Also, we need a simple version of the uniqueness for the system of linear Volterra integral 

equations of the second kind. 

Lemma 2.3. Let n (^ N, a > 0, p > 1, q > 1, 1/p + 1/q = 1. Assume that r 1-^ Wj{t,T) is 
measurable in (0, t) for all t G (0, a] and 



\Wj{t,T)fdT < const., Vt G {0,cr],j = l,n. 



Then the system 



^i(i) = E ilw,{t,T)^,{T)dT j ,Vt G (0,a],i = 1, 



3=^ \0 

has a unique solution {^jYj^i = in L^(0,cr). 

Proof. Using Holder inequality one has 

t 

\Mt)\ < E 



Wj{t,T)'^j{T)dT 



1/p 



1/9 



< 



Y^Umit^Trdr] il\^,{T)\UT 



< con.st.Y^ / \^j{T)\'^dT 



1/9 



Therefore, 



f;i*,(r)r<con.t.f;[ /|vl/,(r)|''<ir| = con.t. /" j f] l^.(r)|M dr, 
j=i i=i \o / \J=i / 

n 

and it follows from Gronwall's Lemma that ^ |^,(t)|'^ = 0. Thus {^j}^^i = 0. D 

Remark 2.3. In later application, we have 

iT^r / M const. . - — 
\Wjit,T)\ <^^, J = l,n. 

\/t — T 

In this case, we can choose any p G (1, 2) in order to apply Lemma W. 



2.2 Reduction to integral equations 

Denote by Gi, G2, G3 the function G in Lemma [2T] corresponding to k = Z), Do,Dj, 

respectively. 

Applying Lemma |2. II to equation (jl.7p with hQ{t) > x > —00 and using condition 

5o(t)/io(i) +DSx{hoit)-,t) = 0, one has 

t /io(0) 

5(a;, t) = -D j Gi^{x, t; /io(r), r)5o(r)dr + ^ Gi(x, t; ^, 0)S{^, 0)d(. (2.3) 

-00 

Thus S{x,t) is determined completely by Hq and So{t). Taking x — > ho{t)— in (12. 3p and 
using the jump relation in Lemma 12.21 we get 

i fto(o) 

So{t) = -2D f Gi^{ho{t),t;ho{r),T)SoiT)dT + 2 j Gi{ho{t),t;tO)S{tO)dC. (2.4) 
-00 

Next, apply Lemma |2. II to equation (jl.6p , for /lo(^) > a: > — cxd, 

t 
;, t) = /" G2(x, t; /io(t), r) [Do«i(t) + To(T)/i^(r)] dr 

' * .0(0) (2.5) 

Do j G2^{x,t;ho{T),T)To{T)dT+ j G2(x,t;e,0)r0(e)de, 



with 

t;i(t)=T,.(/io(t)-,t). 

We differentiate both sides of (12.51) with respect to x, then take x — > /io(i) — • To go into 



the details, because DoG2^x = —DoG2^^ = G2t, we have 

t t 

Do / G2^xix,t;ho{T),T)To{T)dT = - G2T{x,t;ho{T),T)To{T)d7 



d_ 

dr 



{G2ix, t; /io(r), r)) - G2^{x, t; /io(r), r)/io(T) 



To{T)dT 



G2{x,t;ho{0),0)To{0)+ / G2(x, t; /io(r),r)r^(r)dT 



G2x{x, t; /io(r), r)ro(r)/io(r)dr. 

Moreover, 

fto(O) /io(0) 

I G2x{x,t;^,0)T\OdC = - I G2^{x,t;^,0)T^iOdC 

— oo — oo 

hoiO] 

= -G2{x,t;ho{0),0)T\ho{0)) + J G2{x,t;^,0)T^{Od^. 

— oo 

Usmg the compatible condition To(0) = T^{ho{0)), we have, for all /io(i) > x > — cxo 



t 
Tx{x,t) = Do G2x{x,t;ho{T),T)vi{T)dT 



t 



t hoiO) 

+ jG2ix,t;ho{T),T)mT)dr+ J G2{x,t;C,0)T^iOdC. 



(2.6) 



Taking x — > /io(t)— and using Lemma |2. II for the first term, we have 



vi{t) = 2Do / G2x{ho{t),t;ho{T),T)vi{T)dT 



t ho{0) 



(2.7) 



+2 J G2{ho{t),t;ho{T),T)T^{T)dT + 2 J G2{ho{t),t;^,0)T^{Od^. 



Now we consider the heat distribution in false-bottoms. Apply Lemma l2. II to equation 



T3]) . for hu{t) > X > ho{t), 



T{x,t) = Di Gz{x,t;hu{T),T)vi{T)dT 



Gs{x,t;ho{T),T) [DjV2{r)+To{T)h'o{T)] dr 



(2. 



hu{lS) 



DjjGst:{x,t;ho{T),T)To{T)dT+ J Gs{x,t;tO)TyOd(, 



ho{0) 



where 



1.. 



V2{t) = T^{ho{t) + ,t), V3{t) = T^{hu{t)-,t) = ^K{t). 

A/ 
Let us differentiate both sides of ()2.8p with respect to x. We have the same calculation to 

t t 

Di I G3^^{x,t;ho{T),T)To{T)dT = / G3^(a;,t;/io(r),r)To(r)dr 

To{T)dT 
t 

G3{x,t;ho{0),0)ToiO)- /G3(x,t; /io(T),r)r^(r)(ir 



d 

— (G3(x,t;/io(T),T)) - G3^(x,t;/io(T),r)/io(T) 



+ / G3^(x,t;/io(T),r)To(r)/io(r)dr, 



and 
G3.(x,t;C,0)T0(OdC = - I G3g(x,t;C,0)T0(Ode 

ho(0) hoiO) 

hu{0) 

= Gsix,t;hoiO),0)T\hom+ J G3ix,t;^,0)T^{^)d^. 

ho{0) 

Here we have used the compatible condition T^{hu{0)) = 0. Using again the compatible 

condition r°(/io(0)) = ro(0), we find that 

t t 

T^{x,t) = Dj G3:r{x,t;hu{T),T)v3{T)dT - Di / G3^(x, t; /io(r), r)v2(r)dr 



t hu{0) 

G3(x,t;/io(T),r)T^(T)dT+ J G3{x,t;^,0)T^{OdC. 

ho{0) 



(2.9) 



10 



for all hu{t) > X > ho(t). Taking x — > /io(i)+ and x -^ hu{t)— both sides of (j2.9p and 

using (Lemma 12. 2p one has 



V2{t)=2Di I G3:,{ho{t),t;KiT),T)v3{T)dT -2Di I G3^{ho{t),t;ho{T),T)v2{T)dT 

k.m (2-10) 

-2 1 G3{ho{t),t;ho{T),T)mT)dT + 2 I G3(/io(t),t;e,o)r|(eK, 

hoiO) 

and 

t t 

V3{t) = 2Di f G3^{hu{t),t;K{T),T)v3{T)dT -2Di fG3Ahu{t),t;ho{T),T)v2{T)dT 



t h.u{0) 

-2 j G3{hu{t),t-M{T),T)U{T)dT + 2 j G3{hu{t),t;^,0)T^{OdC. 

ho{0) 

(2.11) 
Four equations (12. 4p , (12. 7|) , (I2.10p and ()2.1ip form a system of nonlinear integral equations 
with respect to So{t), vi{t), V2{t) and V3{t). Here, due to the interface condition (jl.ip and 

(|1.4p . one has 

t t 

/io(i) = /io(0) + Xi I V2{T)dT -Xo I vi{T)dT. (2.12) 



and 

t 
K{t) = hu{0) +\lj V3{T)dT. (2.13) 



We thus have proved the direct part of the following statement. 

Theorem 2.2. The problem ( [i.i| )- |777| ) in [0,cr] is equivalent to the problem of finding 
a continuous solution v = {Sq^vi^V2-,V2) in [0,cr] for the system {2.4\) , ( f^.7| ), \2. 1 0\) and 
\2.11\) . where ho and hu are given by ^2.1S\) and \2.13\) and satisfy hu{t) > h^it) in [0,(t). 

Of course, we only need to prove the converse part of the statement. 

Proof. Suppose that v = {So,vi,V2,V2) is a continuous solution in [0, a] for the system 
dMD, (EZD, (I230|) and (l2TT]) . with /iq and /i„ are given by (l2J2]) and (12^31) and satisfy 
Kit) > ho{t) in [0,(t]. Define S{x,t), T{x,t) by ([23]), ([23|) and jSS]). We want to check 
{ho,hu,T,S) satisfies (C1)-(C4). 

It is clear that (CI) is automatically satisfied. Because S{x,t) and T{x,t) are defined by 
integral forms ()2.3p . ()2.5p and ()2.8p . the smoothness conditions in (C2)-(C3), excepting the 
continuity of T{.,t) at x = hQ^t), hold. The initial conditions in (C4) simply follow by 
getting t ^ 0+ in (l2T^ . (I2T3D . (lOD . (l23]) and dZSD and using 

^Um G(x,t;e,0)=(^(x-e), 
11 



with G = Gi, G2, G3. It remains to check equations (|l.ip - (|1.7|) and prove that T{ho{t),t) = 

To{t). 

Step 1. We start by prove three diffusion equations p.5p - (jl.7p . To prove (|1.5p . where 

T{x,t) is defined by (|2.8p in hu{t) > x > /io(i), we shall verify that each of four terms in 

the right-hand side of (j2.8p is a homogeneous solution of the operator (d/dt — Djd'^/dx'^). 

This fact is clear for the fourth term due to the property of Green function, 

(' hu{0) \ 

I Gsix,t;C,0)T''{OdC 
ho{0) J 

hu{0) 

hoiO) 

For the three first terms, we have to be more careful because in general differentiating with 

t 
respect to t a function of the form t >-^ J K(t, T)dT may cause a jump, 





I K{t,T)dT\ = lim K{t,T)+ f Kt{t,T)dT. 



However, in this case the jump lim K(t, r) vanishes due to 

T^t — 

lim G3(x,t;^,r)= lim ^3^(2;, t; ^,r) = 0, Vx / ^. 

Therefore, the three first terms can be treat similarly to the fourth term. Thus (jl.Sh holds. 

And (|1.6p - (jl.7p can be treated by the same way. 

Step 2. Next, we check that 5 fits the Dirichlet condition S{ho{t) — ,t) = So{t), T fits three 

Neumann conditions Tx{ho{t) — ,t) = vi{t), Tx{ho{t)+,t) = V2{t) and Tx{hu{t)—,t) = V3{t), 

and deduce the Stefan conditions (jl.ip and (jl.4p . 

Infact, from (|2.3p taking x -^ /io(t)— and using Lemma 12.21 one has 

t 
Siho{t)-,t) = ^So{t) - D I G,^iho{t),t;ho{T),T)SoiT)dT 



00 



+ y Gi(/io(t),t;e,o)5(e,o)de. 



Comparing the latter equation to (j2.4p . we conclude that S{ho{t) — ,t) = So{t). 
Differentiate both sides of ()1.5p with respect to x and then take x — > hu{t) — . Using the 
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same process of geting (|2.1ip from (j2.8p , one has 

t 

T^iK{t)-,t) = ^V3it) + Dj f G3Ahu{t),t;K{T),T)v3{T)dr 



t 
-Di / G3r,{hu{t),t;ho{T),T)v2{T)dT 


t hu{0) 

- J G3{hu{t),t;ho{T),T)T^{T)dT + I G3{hu{t),t;^,0)T^{Od^. 

ho{0) 

Comparing the latter equation to (j2.11l) . we find that Tx{hu{t) — ,t) = v^{t). Similarly, 

Tx{ho{t)+,t) = V2{t) and T^(/io(t)-, t) = vi{t). Thus the Stefan conditions ^^ and ([lJ]) 

follow the definitions of /io(0 and hu{t) in (j2.12p and (j2.13p . 

Step 3. Finally, we show that S fits the Stefan condition (jl.2p . T fits the Dirichlet conditions 

([OD andr(/io(t),t) =To. 

Applying Lemma [2TT] to equation ()1.7p . for /io(t) > x > — cx) and using S{ho{t) — ,t) = So{t), 

one has 

t 

S{x,t) = j Gi{x,t-ho{r),T) [Z)5,.(/io(t)-,t) + 5o(r)/i^(r)] dr 



t oo 

-D I Gi^{x,t;ho{T),T)So{T)dT + I Gi{x,t;^,0)S{^,0)d^. 

-oo 

Comparing to the original definition of S in (j2.3p . we deduce 
t 

/ Gi(x,t;/io(r),r)^i(r)dr = 0, /io(i) > x > -oo, 


where ^i(t) = DSx{hQ(t) — ,t) + S'o(t)/iQ(t). We next differentiate the latter equation with 
respect to x, then take x — > ho{t)— and use Lemma [2^2] to get 

t 
^i(t) = -2 f Gixihoit),t;ho{T),T)^i{T)dT. (2.14) 



This is a linear Volterra integral equation of the second kind with respect to ^i(i), and 

\Gix{ho{t),t;ho{T),T)\ < '"''^''^' 



It follows from Lemma [O] that ^i(t) = 0. Thus P^ holds. 

Similarly, applying Lemma [2TT] to equation (jl.6p . for /io(i) > x > — oo, one has 

t 

T(x,t)= /"G2(x,t;/io(r),r)[Z)or,(/io(r)-,r)+r(/io(r)-,r)/i'o(r)]dT 



t ho{0) 

- Do I G2i^{x,t;ho{T),T)T{ho{T)-,T)dT + j G2{x,t-i,Q)T\i)di, 
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Comparing the latter equation to the original definition of T in (j2.5p . and using Tx{hQ{t)—, t) = 
vi{t), we obtain 

t 
J [G2{x, t; ho{T),T)h'o{T) - DoG2^{x, t; /io(r), r)] ^2{r)dT = 0, (2.15) 



where ^2{t) = T{ho{t)-,t) - To{t). Taking x -^ K{t)- in ([213]) and using Lemma [221 we 
find that 

t 

^2(t) = 2 / [G2{x,t;ho{T),T)h'o{T) - DoG2dx,t;ho{T),T)] ^2{T)dT. (2.16) 



This is a linear Voterra integral equation of the second kind. Since 

const. 



|G2(x,t;/io(r),r)/io(r) - DoG2i^{x,t;ho{T),T)\ < 



Vt 



it follows from Lemma [2.3l that (|2.16p has a unique solution ^2(i) = 0. Thus T{hQ{t)—, t) = 

To{t). 

Applying Lemma [2. II to equation (|1.5p . for hu(t) > x > hQ{t) and comparing to the original 

definition of T in (|2.8p , one has 



f [G3{x,t;ho{T),T)h'QiT) - DjG3^{x,t;ho{T),T)]^3{T)dT 



(2.17) 
[Gsix, t; hu{T),T)h'^{T) - DjG^^ix, t; /i„(r), r)] *4(T)dT, 

where ^3(t) = T{ho{t)+,t) - To{t) and ^4(t) = T{hu{t),t). Taking x ^ /io(i)+ and 
a; ^ hu{t)— in (|2.17p . we have 

t 

[t] = J [G3{ho{t),t;ho{T),T)h'o{T) - DjG3dho{t),t;ho{T),T)]^s{r)dT 





t 
- j [G3iho{t),t;huir),T)h'^iT) - DiG3^ihoit),t;huiT),T)]^4{T)dT, 



and 

t 

^4(t) = - / [G3{Kit),t;ho{T),T)h'o{T) - DjG3i{hu{t),t;ho{T),T)] *3(r)dT 



(2.18) 



(2.19) 



+ / [G3{K{t),t;huiT),T)h'^iT) - DiG3^{Kit),t;hu{T),T)] ^4(r)d7 




Equations (|2.18p and ([2.19P form a system of linear Volterra integral equations of the second 
kind, and we deduce from Lemma 12.31 that ^3 = ^4 = 0. Thus ([l.Sp holds, and the 
continuity of T{., t) at ho{t) follows T{ho{t)+, t) = T{ho(t)-, t) = To{t). D 
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Remark 2.4. The technique in Step 3 is learnt from A.Friedman J2j. However, there is a 
minor mistake in 12]: formula (1.28) in page 221, i.e. 



Gc{x, t; s{t),t)u(s{t), T)dT = 0, 



should be 



[G^{x, t; s{t),t) - G{x, t; .s{t),t).s' {t)] u{s{T),T)dT = 0. 



Fortunately, the conclusion u{s{t),t) = still holds by the same argument. 

2.3 Solving integral equations 

Now we accomplish the proof of Theorem 12.11 We want to prove the system (|2.3p , 
()2.7p . (|2.10p . (j2.1ip has a unique local solution, and this solution can be extended uniquely 
by prolongation whenever the condition hu{t) > ho{t) still holds. The system (12. 31) . ()2.7p . 
()2.10p . (j2.1ip can be rewritten as v = Pv, where 

V = {So,Vi,V2,V3),Pv = {PiV,P2V,P3V,PiV) . 

Denote by C{a,M) the space of functions v = {So,vi,V2,V3) continuous in [0,cr] and 
IbllfOal := sup max{\So{t)\,\vi{t)\,\v2{t)\,\v3{t)\} < M. 

t6[0,(7] 

We shall prove that for M > large enough, says M > 2 ||-P(0)||rocri5 there exists tM > 
depending on M (but independent on w S C{tM,M) ) such that P is a contraction on 

C{tM,M). 
First at all, note that for each M > 0, if t^ > small enough then it follows from 

(l2T2D - (l2Ti|) that 

/i„(0) - /io(0) 



hu{ti) - ho{t2) > 



>0, Vti,t2G [0,tM]- 



(2.20) 



In particular, hu{t) > hQ{t) for all t E [0, ^^1) and hence P is well-defined. 
Let us estimate \\Pv — Pi^llrQ^ i for v,v £ C{tM,M). In what follows, denote by h^ and 

ho the functions given by (J2.12p - (j2.13p where v is replaced by v; and denote by Co > a 

constant depending only on data (/io(0), /im(0),T'^(x), 5'^(x),To(t)) but independent on M, 

V and V. We shall go to the details for P4, and Pi,P2, P3 can be treated by the same way. 

For preparation, we need some estimates related to Green's function. 

Lemma 2.4. Let v,v £ C{tM,M). Then for t > small enough we have 



G3x{huit),t; hu{T),T) - GsxiKit),!; hu{T),T) 

G3{hu{t),t-ho{T),T)-G3iK{t),t;ho{T),T 



GoM „ 



< t\\v — v\ 



G3xihu{t),t;ho{T),T) - G3^(/i„(i),t;/io(r),r) 



I [0,tM 

< t\\v — V 



[0,tM] 



G3xihu{t),t;^,T) - G3^{hu{t),t;i,T 



< 



<Go 



exp 



4Dit 



+ exp 



ADit 



V — V 



[O.tjv 



for all t G [0, ^m] and ^ eR. 
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Proof. First at all, we note that /iq ,/iu are linear with respect to v and 



\hu{t) - hu{T)\ 



Xi / V3{r)dt 



<^l\\v3\\[0,tj,,]it-T), 



|/io(t) - /io(t)| < [Ao ll^i||[o,t„] + Ao ||^2||[o,t„]J (t-r). 
Let us now consider function 

G3x{K{t),t;hu{T),T) = ; =exp 



as a one- variable function of {hu{t) — huir)). Using Lagrange formula one has 

G3x{hu{t),t;huiT),T) - G3a:ihuit),t;hu{T),T) 



{K{t) - K{t)) - [K{t)-K{T) 



exp 



2^2 



The first inequality follows \9\ < XjAI{t — r) and 

{hu{t) - huir)) - \hu{t) -huir)] < \i{t - t) \\v -vW^q^Im] 
By the same argument, one has 

G3(/i„(t), t; /io(r), r) - G3{hu{t),t- ho{T),T) 

K{t)-hu{t)\e 

■ exp 



4Di{t - t] 



AjnDnt-rf 'V ^Djit-r)J^ 



and 



G3x{hu(t),t;ho{T),T) -G3x{hu{t),t;ho{T),T) 
hu{t) - K{t) 



: exp 



ADl{t-T] 



29^ 



4^7rL'f(t-r)3 

for some 9 between {hu{t) — holr)) and {hu{t) — ho{T)). Note that 

\huit) - hiu{t)\ < XiM \\v - wi||[o,i] , 
and it follows from ()2.20p that 

/i„(0) - /io(0) 



ADi{t - t] 



> 



>0. 



Therefore, the second and the third inequalities follow the elementary inequalities z^e~^ < 
const., and {z^ + z'*)e^^ < const., for z > 0. 
For the last inequality of Lemma 12.41 we write 

G3{hu{t),t;C,0) - G3{hu{t),t;C,0) 

exp — ^ — exp ' 



2^/TTlJrt 



ADi{t-T) 
16 



ADi{t-T) 



exp(-|]+exp(-|- 



and then employ the following elementary inequality 

|exp(— Z;^) — exp(— ^2)! < const 

for zi,Z2 G M. To prove (j2.2ip . we can assume < zi < Z2 and note that 

exp{-zj) - exp(-z|^ 



Zl - Z2\ 



(2.21) 



Zl - Z2 



-20exp(-r 



for some 6 G {zi,Z2- Then p.2ip follows |0| exp(— 6*^/2) < const., and exp(- 
exp(-z|/2). 



?72) < 
D 

Return to estimate \\PiV — P4'y| | [o,t^,j] • Using the first inequality of Lemma \2M we get 
G3x{huit),t; huiT),T)v3{T) - G-ix{hu{t),t] /iu(r), r)w3(T) 

G3x{hu{t),t; hu{T),T) - G3x{hu{t),t; hu{T),T) \v3{t)\ 
+ \G3x{huit),t;hu{T),T)\ \v3{t) -V3{t)\ 



< 



^ CqM „ ,,, , CqM 



Vt 



Therefore 



21)/ 



G3x{hu{t),t]hu{T),T)v3{T) - G3xiKit),t; K{t),t)v3{t) 



dr 



Co(M^ + M) „ 
< 2^)/ / ^ \\v — v\ 

^t-T 



\Q,tM 



dr = ADiCo{M^ + M)./t^ \\v - v\\ tq , , . 



Also, from the second and the third inequality of Lemma l2.4^ we have 



2Di 



G3x{hu{t),t; /io(t), r)f2(T) - G3x{hu{t),t; /io(r), r)w2(T; 

t 
< 2Di j t{M + 1) \\v - v^,t,,^dT = 2Di{M + 1)4 ||^ - v\\^^^^^ 



dr 



and 



2Di j \G3{K{t),t; ho{T),T) - G3{huit),t- /io(r), r) | r^(r)d7 


t 

<2Di / t\\v-v\\^Q^t^^]T^{T)dT<Cot\\v-v\\^Q^^^^y 



(2.22) 



(2.23) 



(2.24) 
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It follows from the last inequality of Lemma 12.41 that 



huiO) 



ho{0) 






< Cn \\v — v\ 



[0,ti( 



/io(0) 



exp 



{i-K{t)y 

Cat 



+ exp 



Cot 



dt 



On the other hand, by making some change of variables, one has 



exp 



Cot 



+ exp 



{^-hu{t)y 

Cot 



di 



oo 

2 J exp(-^\d^ = 2y^^^t\ 



V — VWlr.. 



Thus 



hu(0) 



C3iKit),t- e, 0) - G3iKit),t; e, 0) T^iC)dC 



ho{0) 



(2.25) 



<2A/^C|iM|k-t'||[o,t^,] 



It follows from (l2:22]l - (12:25] ) that 

\PaV - PAv\\\ot..] < CoM'^^/t^\\v -v\\ 



I [0,tAf ] 



[0,t], 



where Co stands for a constant depending only on data {ho{0),hu{0),T^{x), S^{x),To{t)). 
We also have similar estimates for Pi,P2,P3. Thus, if 

M^VU^ < eo, (2.26) 

where eo stands for a constant depending only on data (/io(0), /iu(0),T*^(a;), 5*^(x),To(i)), 

then 



1 



^^-^^'lllo.tM] ^ 9 ll^-^ll[o,tM]' Vv,f GC(tM,M). 



(2.27) 



The latter inequality implies that if M > 2 ||P(0)|| fq i then P maps from C{tM, M) into 
itself and is a contraction. Thus P has a unique fixed point v in C{tM, M) for tM small 

enough (restricted only by (|2.26p ). 
Now let us prove the solution can be extended uniquely by prolongation. Assume that v is 
a solutions of system 1^, ([221), (l2lnD . (f2TT]) in [0, a*] with a* < a and K{a*) > ho{a*). 
Let us consider the system ([23]), (122]), fl2TIiD . (1211]) for t > a* (instead of t > under the 

data 

{ho{a*), K{a*),T{x, a*), S{x, a*),To{t)). 

We need to check that the new data satisfy Hypothesis (H1)-(H4). Infact, (HI) and (H4) 
are automatically satisfied because hu{cr*) > ho{cr*) and To(t) is even continuously 
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differentiable in [0, a]. For (H2), the continuity of T(.,cr*) and Tx{-,a*) is guaranteed by 
the definition of the solution. The boundedness of Tx{-,(t*) fohows the finite hmits 

T^{ho{a*)+,a*)=V2{a*), 

T,{hu{a*)-,a*) = V3{a*), 

and, from (|2.6|) . 

ho{0) 

G2(x,t;e,o)r|(eK 



hm sup 



hmsup \Tx{x, a* 

Here we have used a property of Green's function 



< sup |tO(0| 



J \G2{x,t-^,0)\d^ = l. 



Similarly, (H3) holds. Thus Hypothesis (H1)-(H4) holds for the new data 
{ho{a*),hu{a*),T{x,a*),Six,a*),To{t)). 

Therefore, it follows from the above result that the solution can be extended uniquely in 

[0,(T**] for some a** € {a*, a]. 

Remark 2.5. It is sometimes difficult to m.easure the initial condition (T(x,0), 5(x,0)) in 
the infinite domain /iu(0) > x > — oo as in the previous section. However, this condition is 
maybe replaced by the initial condition (r(x,0), 5(x,0)) in a finite interval hu{0) > x > L 
and the history information at one point {T{L,t), S{L,t)) inO <t < a, where L is a deep 
point in the ocean satisfying L < /io(i) for < t < a. 
Infact, introduce the Green's function for the half-plane x > L, 

Kj{x,t;(,T) = Gj{x,t;C,T) - Gj{2L - x,t;C,T), j = 1,2. 

Note that Kj{x,t;L,T) = and the formulae in Lemma \2.1\ and Lemma \2.2\ still hold with 
Kj replacing Gj and a{t) = L. Thus instead of Ii2.cl\) and Ii2. 5\) we have, for /io(i) > x > L, 
t t 

S{x,t) = -D f Ki^{x,t;ho{T),T)So{T)dT + D f Ki^{x,t;L,T)S{L,T)dT 

° (2.28) 




ho(0) 



+ J K,ix,t;C,0)S{C,0)dC, 

L 

t 

Tix,t) = J K2{x,t-ho{T),T)[Domho{r)-,T)+ToiT)h'oiT)]dT 



t t 

-Do I K2^ix,t;hoir),T)To{T)dT + Do I K2^ix,t;L,T)T{L,T)di 



hoiO) 

+ I K2{x,t;^,0)T{(,0)dt 



(2.29) 
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Therefore, the uniqueness of the solution {hQ,hu,T, S) of the system il.l\) - (T77\ ) can be 
proved by the same argument. 
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